METHOD OF ANALYZING ELECTROMAGNETIC INTERFERENCE 

BACKGROUND OF THE INVENTION 
l,Field of the Invention 
5 The present invention relates to a method of analyzing electromagnetic 

interference (EMI) (hereinafter often referred to as an "EMI analysis method"), 
and more particularly, to a method of analyzing electromagnetic interference 
arising in a large-scale, high-speed LSI (large-scale integrated circuit) by means 
of high-speed, highly-accurate logic simulation* 
10 l.Description of the related Art 

LSIs fad — are used in a broadening range of applications, from 
communications devices, such as cellular phones, to general household products, 
toys, and automobiles, as well as applications in the field of computers. 
Electromagnetic interference arising in such — a products induces radio 
15 interference noise in a nearby RF receiver, such as a TV set or a radio, or in 
faulty operations of anothe r nearby system. In order to prevent this problem, 
the entirety of a whole p roduct is shielded, or filters are provided in a — the 
product. With a view io^^sf-diS order to reduce tbe need to 

increase iB-the number of comoomnts comi>Iexity a nd cost of oroduction a nd to 
20 reduce the difHeulty difficu^^^^ encountered i n coiiventfonal rnethods of 
preventing occurrence of electromagnetic interference in a product, strong 
demand exists for a means of suppression of noise in an LSI itself. 

linder SH«4> a i5lt^iatioj>j sAn LSI is f-af^ked-as a key flevi<^<^ f oi-com^mrent 
of any product in which €oaf.ai ff >S"£in LS l it is used . Demand exists for a^^ 
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larger-scale, high-speed LSI for cyis^rmg l n order to eiisftii'e the competitiveness 
of a product. In a s itu atioii in which the cycle oFA s product development cycles 
l>eeomes- - become s horter, design- automation of aiv -LSIs is— has bccoiii e 
indispensable for satisfying the demand. There is a growing necessity n eed for 
5 adopting synchronous circuit ry in LSIs as a coadMoR for mtroduci ng a i n order 
to take advantage of state-of-the-art design-automation methods. In a case 
where When a ll of the circuits of a large-scale, high-speed LSI operate 
synchronously with a reference clock signal, instantaneously-changing current 
becomes very large and induccsi resulling in an increase in the indu ction of 

10 electromagnetic interference. 

The present invention relates to a simulation method which that e nables 
evaluation of EMI indispensable for reducing electromagnetic interference while 
maintaining a tendency toward a larger-scale, higher- speed LSI. 

Noise imposed on another device by an LSI is roughly classified into two 

15 typesti radiation noise, and conduction noise. Radiation noise emanated 
directly from an LSI includes noise emitted from internal wires of an LSI. 
However, the i nternal wires di> -are n ot act as an antenna of l arge sigc enough to 
act as an effective antenna * As oneratiiig frequencies of LSIs continue to 
increase, a matter of course, i t is cons i dered expected t hat the radiation n oise 

20 emitted directly from an LSI will pose a problem in the futur e^ in iissoeiatlen 
with an improvemeot in the operation frequency of an LSI . However, as of 
now, the noise emitted from the inside of an LSI is considered trivial. 
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In contrast, conduction noise from an LSI a ffects smother devices 
mounted on a -the same p rinted circuit b oard, by way of direct interconnections, 
such as w4r-€^" <?\t€mai leads of af> "the L SI or i-oii^ln^ -traccs p rovided on a- tbe 
printed v¥irin$ £ circuit board> Noise is emitted from sucii interconnections 
5 while t The interconnections a^^e-act as the source of origination or as an antenna 
for emittiug conductioa noise . The antenna constituted of the interconnections 
is much larger than that constituted by internal wires of an LSI and is a 
dominant element in terms of electromagnetic emission. 

A power line and a signal line can act as paths along which conduction 

10 noise developing in an LSI travels. In consi der ati o n o F When analyzing an 
electromagnetic field in the vicinity of an LSI, noise which results from 
variation in an electric current of a power source being emitted from a power 
line serving as an antenna is considered to be dominant. There may 1^ a In 
some c ases, where ringing and overshoot phenomena stemming from variation 

15 in a signal may also p ose problems. However, theye more frc<|ucntly arlse& a 
case where variation in an internal power level of an LSI J^at»propagates as a 
signal waveform, to ther e by more often p resents a problem. Noise emitted 
from a power line or a signal line is considered to have a strong correlation with 
variation in the electric current of a power source (hereinafter referred to as a 

20 "source current'')* 

A source current of a CMOS circuit will now be described by reference 
to a simple inverter circuit. In a case where When variation arises in a voltage 
applied to an inverter circuit, ther^- flow^ a load capacity charge/discharge 
current flows, which is the primary source current of the CMOS circuit. In 
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addition, a short circuit current flows together with the load capacity 
charge/discharge current. In automated d esign of such a CMOS circuit, all 
circuits of an LSI are synchronized in accordance with the c onstraints on: the 
esc of a— liie design-automation tool. As a result of all circuits being 
5 synchronized, all circuits of the LSI operate simultaneously, and a peak current 
arises in a power source in synchronism with a reference clock signal. Further, 
in order to increase operating speed, or shorten a cycle, of the LSI, the capacity 
of a transistor is increased so as to enable a charging/discharging operation to 
be completed within a shorter period of time. Eventually, a-the peak current 

10 increases. As^-a matter of coursc N ecessarilY^ the total source current of an LSI 
is increased when the integration l evel of an LSI is increased. Thus, the peak 
current of the power source is increased, thereby inducing occurrence of a n 
abrupt change in a source current. Such an abrupt change induces an increase 
in higher harmonic components, thereby resulting in an increase in 

1 5 electromagnetic interference. 

Ili^hl v - i>rceis e Highl v precise simulation of change in a source current, 
which may be said to t he p rimarily account for cmse of e lectromagnetic 
interference, is considered to be effective in ev-alHi$:tio3^--^f - predicting the 
electromagnetic interference arising that will arise i n an LSI. 

20 A conyentional c urrent simulation method for eiTeeting transistor-level 

current analysisr-as will be described below , has coovcntiojiia il y been employed . 

FIG. 15 is a block diagram showing the flow of processing operations 

peFt-armif^ to of a conventional transistor-level EMI analysis method. 

According to this method, en - fhe fo^isi sr o f based on layout inforgnarien data is 
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provided 01 that D crtamiiie to describes a n LSI ^vbieh that i s to be analyzed 

th ro ugh use of t isins a transistor-level current analysis method, ^ there is 

i>erformed 4a¥Out Lav ou t p arameter extraction (hereinafter referred to simply 
as an "LPE'O processing 03 is performed on the layout data . Subsequently, 
5 there is several processing steps are p erformed: circuit simulation 06 regarding 
of a switch-level netlist; source-of-current modeling 08; a power line LPE step 
OlO; transient analysis simulation 012; and fast Fourier transformation 
(hereinafter abbreviated F FT) processing 014. 

Processing pertaining to each of the foregoing processing steps vdll now 

10 be described by"W|lk,reference to FIG. 15. 

First, in step 03 data a ^e-is input: layout data Ol pertaining to a 
semiconductor integrated circuit to be subjected to EMI analysis; parameters of 
elements, such as transistor elements or various parasitic wiring elements (e.g., 
resistors and capacitors); and an LPE rule 02 for defming a form in which 

15 extracted layout parameters are to be output. In accordance with the LPE rule 
02, parameters of the respective elements included in the layout data Ol are 
calculated, whereby a netlist 04 is produced. In step 03, parasitic elements of 
a power source (and the ground) are not ofeieets-^f eMt?;=aetloi^ ext?'acted . 

In step 06 arc input t he netlist 04 prepared in step 03 and a test 

20 pattern 05 are ia put. The test patt era is used for causing a -the c ircuit bein^ 
analyzed , which ser%^es iis an object of aimlvsis. to replicate a desired logic 
operation. There are calcuifttcd a A load capacity charge/discharge current 
and a short circuit current arc catcuiatcd. which correspond to the operating 
state of an internal circuit, thereby producing current waveform information 07 
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€OH€€rHmg - a bout the waveform of an electric current of a transistor* The first 
op ^ ei - ati^R of th e p rocessing pertaining to of-Step 06 is e ffected based o n the 
assumption that the potential of a power source (and that of ground) is a 
variation- free, ideal potential. 
5 In step 08 is ent < ei -e d t he current waveform information 07 eoaccrnii^g 

of a transistor prepared in step 06 is entei-ed . The thus entered current 
waveform information 07 is modeled into a mode which can be applied to 
subseeueHtH>tei> -Ot2v--w-hef^wltj^^^ ^ to prepare c urrent source element model 
information 09 is prepared saitafoie for subseiiuent sleg 012 . In order to 

10 alleviat e a reduce the processing load which would be imposed on f or 
subsequent step 012, a function circirit block consisting of a plurality of 
transistors is usually modeled as a single current-source element. 

The F i'ocessing- processing i> ert-a>Hlf>g-to- i?erform€d in s tep OlO differs 
from processing pertaining to -step 03, only in that rather than parameters of 

1 5 transistor - elements - and those - of - varlous - p 

of parasitic elements of a power source and those of a ground wire (e.g., 
resistors, decoupling capacitance, and like elements) ar c tal^cn as objects of 
extraetlo ^ extracted^ rather than parameters of transistor elements and of 
various parasitic wiring elements . Hence, repeated explanation is omitted. In 

20 step OlO, a power source (and ground) wiring netlist Oil is produced. 

In step 012 ore eiiter^- the current source element model information 
09 prepared in step 08„ijS„ entered, the power source (and ground) wiring netlist 
Oil prepared in step OlO is entered, and impedance 016 of a wire or a lead 
frame (including, resistance, capacitance, and inductance) is entered . Through 
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analysis of these input data carried out by a transient analysis simulator typified 
by SPICE, fluctuations in line voltage of a circ uit to be analyzed subject circuit 
are calculated. Thus, tfe<H-e 4^ pi^duccfi a line voltage drop result 017 
concerning — the — is i^roduced wh ich co r resiionds to thus - calcuiQted — these 
5 fluctuations in line voltage. 

Subsequently, the p rocessing oertaimng to of s tep 06 is performed 
again. In contrast with the first operation of the processing pertaining to o f 
step 06 hawng""been---effeeted-- -w was based o n the assumption that the 
potential of the power source (and the ground) is a fluctuation-free, ideal 
10 potential, the line voltage drop result 017 prepared in step 012 is entered. 
The current waveform information 07 coiiceriiing f or a transistor is prepared 
again m — with c onsideration of fluctuations in Une voltage. Similarly, 
processing peytai?>mg -to^-of steps 08 and 012 is repeated. 

Processing pcrtuining to steps Steps 06, 08, and 012 is effected a re 
15 repeated several times in a l oo p e d manner, wherewith there - is produced thei^b^^ ^^ 
producing a current waveform result 013 wliieh that very Mgbly - a ccurately 
duplicates simulates fl uctuations in line voltage. 

In step 014, the current waveform result 013 prepared in step 012 is 
entered and subjected to fast Fourier tyansforii^at i on {hereinafter abbrcv i at 
20 FFT) - FFT processins ^ to thereby enable frequency spectrum analysis. - Th e re -i s 
obtained ^ luSa^an EMI analysis result 015 is obtained . 

In the conventional example, the precision of verification varies greatly 
according to the combination of the LPE processing 03, the power line LPE 
processing OlO, and the source current modeling processing 08. However, a 
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certain level of accuracy of analysis can be expected. A transient analysis 
simulator typified by SPICE is Hsed f<>r l imited to t ransistor-level analysis of an 
electric current.- — MeiieejHa4mHtat4mi4s-4mpo 

a nalyz e d, and thus^ an enormous amount of processing time is required. 
5 Since the Integration level of a~^emiconductor integrated circuits has increased 
recently, ftBd--establishment of an EMI analysis method which ' i s desired to 
enablea^ high-speed analysis of an electric current on a higher level larger t han a 
transistor level-is-desi?=ed. 

A gate-level current analysis method has conventionally been proposed 

10 as a current analysis method wMeb t hat c an be made faster. This gate-level 
current analysis method is used for analyzing power consumption. One 
example of a gate-level current analysis method is EMI-noise analysis which is 
%e effeeted- that is oerfoimed i n an ASIC design environment, l-fee -rhis 
method is described in "EMI-Noise Analysis Under ASIC Design Environment" 

15 (ISPD&99, pp. 16 through 21). According to this technique, an event is 
acquired from the result of a gate-level simulation using a test vector, and the 
waveform of an electric current is estimated. The frequency of the 
thus- estimated current waveform is analyzed tferm^Efa "•fast"" -iismg F oi^ri<^r 
Iranstormation (F FT 4 processing . More specifically, as shown in FIG. 16, a 

20 logic simulation 104 is egeeted on th e- b^i^ - of ba^^ netlist 101 and a test 

vector 102, wherewith event information 105 is calculated. On the basis of the 
thus calculated B ased on this e vent information 105 and on waveform 
information 103 obtained at the time of toggling, processing p<H-t4m»ng--tO"£or_a 
current waveform calculation section 107 is executed, to thereby produce 
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nroducing a current waveform calculation result 108. This current waveform 
calculation result 108 is subjected to FFT processing 109y to th e reby p roduce a 
frequency characteristic 110. The EMI-noise analysis method can elf^t 
perform a n EMI analysis operation faster than that performed according to t he 
5 conventional gate — ^level EMI analysis- method. However, use of a test vector 
still involves consumption of much substa alia! e xecution time. Therefore, the 
processing speed achieved by the aforemen tioned EMI-noise analysis method is 
not sufficiently high, and demand still exists for an- iim-ease iiv-p^i-o<^essing"Speed 
of the f aster E MI-noise analysis method. Another problem ^ I^e — with the 
10 aforementioned E MI-noise analysis method js Qiat al i so eneou a t ers a - problem of 
aa -the a nalysis result bem^is dependent on the pattern of an emoioved the test 
vecto r employed . 

As mentioned above, the conventional example using the transistor-level 
current analysis method can be expected to yield a certain level of accuracy. 

15 However, a transient analysis simulator typified by SPICE is used for such a 
transistor-level current analysis. As such. A r^a Umitation is imposed on the level 
of a circuit to be analyzed, and an enormous amount of processing time is 
required. The level of a- semiconductor integrated circuits has recently 
increased recently , and thus t here is desired estublishment of a need for an EMI 

20 analysis method wMch method that enables high-speed analysis of an electric 
current at a scale larger than the scale that which can be analyzed by a 
transistor-level simulator. 

The gate-levi?! G ate- level s imulation using a test vector has also been 
proposed. However, the example conventional gate-level simulation technique 
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encomiters ditficoltv in increa s ms does not sufficiently mcrease t he speed of 
analysis. Since the gate-level simulation technique employs a test vector, an 
analysis result is dependent on the employed test patter n employed * 

SUMMARY OF THE INVENTION 
5 The present invention has been conceived to solve the drawbacks of the 

conventional methods and is aimed at evaluating electromagnetic interference 
developing in an LSI through a simulation by means of high-speedy 
highly-accurate analysis of a power-supply current* 

To this end, the present invention provides a method of analyzing 
10 electromagnetic interference (an EMI analysis method). In contrast with a-the 
known dynamic gate-level simulation method, the present E MI analysis method 
enables estimation of EMI noise^ by means of calculating signal propagation of 
each node ttomi^ -me"<»f -using t he signal propagation probability technique, 
and calculating variation time of each node through use of u sm^ ^^he Static 
15 timing analysis technique^\ In short, the present invention is characterized in 
involves c alculating a frequency characteristic from the relationship between 
toggle probabiUty of each node and delay in each node. 

More specifically, the present invention provides a method of analyzing 
electromagnetic interference generatin g developing in an LSI, comprising: 
20 a correction step of correcting the amplitude of a current estimation 

waveform in each node which has been previously prepared for each change in 
each node, in accordance with the probability of variation in each node; 

an addition step of adding current waveforms of all nodes together 
within a period of time corresponding to one cycle, provided that the 
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thus-corrected current waveform appears at a time a signal arrives at each 
node; and 

a frequency analysis step of analyzing the frequency of the current 
waveform calculated in the addition step. 
5 The probabiUty of variation in each node is calculated th rough use of 

using t he signal propagation probability technique. Further, the time at which 
a signal arrives at each node is calculated through use of u sinj^ the static timing 
analysis technique. The time at which a signal arrives is defined so as to fall 
within the range between the maximum time and the minimum time, in 
10 accordance with an average time, the maximum time, the minimum time, or a 
predetermined distribution such as a normal distribution. 

Preferably, in the correction step, the ampUtude of a current estimation 
waveform, which has been prepared for each change in each node, is corrected 
in accordance with the probability of variation in each node and a distribution 
15 with respect to time (hereinafter called "chronological distribution'^). 

More specifically, the probability of change in each node is calculated 
throu^ use of u sing t he signal propagation probability technique, and the 
chronological distribution at which a signal arrives at each node is calculated 
Ihrou^^ use of usins^ t he static timing analysis technique. 
20 Preferably, each node has a pluraUty of paths, and a current waveform 

is calculated in consideration of based on a case where each of the paths has a 
unique probabiUty of change and signal arrival time. 

Further, the present invention provides a method of analyzing 
electromagnetic interference developing in an LSI, the method comprising: 
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a waveform formation step of forming a current estimation waveform 
which has been prepared for each change in each node^ as if the waveform 
randomly arises within a plurality of predetermined cycles, in accordance with 
the probability of change in each node and a time at which a signal arrives at 
5 each node; 

adding the thus-prepared current estimation waveforms of all nodes, to 
thereby derive a current waveform; and 

analyzing the frequency of the current waveform, thereby determining a 
noise characteristic of EMI. 
10 In other words, under the EMI analysis method according to the 

present invention, the probability of change in each node is calculated through 
use of using t he signal propagation probability technique, and the result of 
calculation is stored as a probability at which a signal randomly changes. 
Further, a time at which a signal arrives at each node is calculated through ^ii^ 
15 frf ^using t he static timing analysis technique. 

Moreover, the present invention provides a method of analyzing 
electromagnetic interference developing in an LSI, the method comprising: 

a waveform formation step of forming a current estimation waveform 
which has been prepared for each change in each node, as if the waveform 
20 randomly arises within a plurality of predetermined cycles, in accordance with 
the probabiUty of change in each node and chronological distribution 
probabiUty; 

adding the thus-prepared current estimation waveforms of all nodes, to 
thereby derive a current waveform; and 
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analyzing the frequency of the current waveform, thereby determining a 
noise characteristic of EMI. 

More specifically, the probability of change in each node is calculated 
rhrou^ us e of u sing the signal propagation probability technique, and the 
5 result of calculation is stored as the probability of a signal changing randomly. 
A chronological distribution at which a signal arrives at each node is calculated 
through use of nsin^ t he static timing analysis technique. 

Preferably, each node has a plurality of paths, and a current waveform 
is calculated in consideration of b ased on a case where each of the paths has a 
10 unique probability of change and signal arrival time. 

BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. 1 is a schematic illustration for describing of the concept of the 
present invention; 

FIG. 2 is a block diagram showing a portion of a circuit used in a first 
15 embodiment of the present invention; 

FIGS. 3A and 3B are graphs showing waveforms of signals arriving at 
the respective nodes of the cell shown in FIG. 2; 

FIG. 4 is a block diagram showing the processing of a frequency 
characteristic calculation block according to a first embodiment of the present 
20 invention; 

FIGS. 5 A to 5D are illustrations for describing s howing a processing 
image according to the first embodiment of the present invention; 

FIG. 6 is a flowchart of processing of a current waveform calcidation 
according to the first embodiment; 
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FIG. 7 is a block diagram showing the processing of a frequency 
characteristic calculation block according to a second embodiment of the 
present invention; 

FIGS. 8A through 8D are illustrations for describing s howing a 
processing image pertaining to the second embodiment; 

FIG. 9 is a flowchart of processing of a current waveform calculation 
according to the second embodiment; 

FIG. 10 is a block diagram showing a portion of a circuit used in a third 
embodiment of the present invention; 

FIG. 11 is a graph showing signal waveforms of each path; 

FIG. 12 is a block diagram showing a frequency characteristic 
calculation block according to a third embodiment of the present invention; 

FIGS. 13A through 13C are illustrations showing a processing 
pertaining image to the third embodiment; 

FIG. 14 is a flowchart of current waveform calculation processing 
according to the third embodiment; 

FIG. 15 is a flowchart for describing s howin<^ a known EMI analysis 
method; and 

FIG. 16 is a flowchart for describing s howing a gat€*level EMI dvnauuc 
analysis method of analyzing EMI dynamic - at gate-leyei . 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

An electromagnetic interference analysis method according to preferred 
embodiments of the present invention will now be described bv- with r eference to 
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the accompanying drawings. As shown in FIG, 1, an EMI analysis method 
according to the present invention i s char actcrigcd i n comprises : 

calculating the transition probabiUty of a node from a netlist 1 and a 
transition probability 2 , Ihroiigh nsc o f ji sing a propagation probability method, 
5 and calculating a static delay 4 throHgh" - iiS<? of u sin^ a static delay analysis 
method, to thereby derive a calculated probability/delay 5 of a -the n ode; 

estimating the waveform 6 of an electric current on the basis of based 
on the probability/delay 5 and information 3 eoH^eynti^g about t he waveform of 
an electric signal at the time of toggling, to thereby derive a current waveform 
10 estimation result 7; and 

subjecting the current waveform estimation result 7 to a fast Fourier 
transformation (FFT) 8 (hereinafter called an ^TFT^Q , thereby determining a 
frequency characteristic 9 of the waveform, 
(First Embodiment) 

15 A method of analyzing electromagnetic interference according to a first 

embodiment of the present invention will be described hereinbelow. As can be 
seen from a The s chematic diagram shown in FIG. 1, imder shows a n EMI 
analysis method according to the pi=ese?>t"fi^t_c™l><>diiiicntv^ 
ibe -The ouantitv of electromagnetic interference developing in an LSI is to be 

20 analyzed on the bag^i ? j •■■■ of b asgd on a transient probabiUty and static delay 
propagation data , provided that a . A waveform shown in FIG. 3A appears ks-at a 
node A of a flip-flop (FF) cell and a waveform shown in FIG. 3B appears m-at a 
node B of the FF cell (wS^ere FIGr aa^nJayged v^ 

the signal designated by braces in FiC. 3A) when a clock signal CLK is input to 
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a circuit shown in FIG. 2 (where FIG. 3B is au eaiarged view of about 1>5 cycles 
of the signal designated by braces in FIG. 3A) . Here, the transition probability 
of a node is calculated from a previously-prepared netlist 1 and a transition 
probability 2. Further, a static delay 4 in a current estimation waveform per 
5 change is calculated. The ampUtude of a current waveform is corrected 
consideration of based on i nformation 3 concermng aboxit t he waveform of an 
electric current arising at the time of a predetermined toggling operation. 
Provided that the corrected current waveform arises at a time at which a signal 
arrives at the respective node, the current waveforms which appear at all nodes 

10 during a period of time corresponding to one cycle are added to the current 
waveform (the current waveform estimation processing 6). The current 
waveform estimation result 7 determined through addition is subjected to the 
FFT processing 8, thereby determining the frequency characteristic 9 of EMI 
components of a circuit to be analyzed. 

15 FIG. 4 is a block diagram fo r des cribing s hoeing the overall prossecing 

flow of proce s sing of the EMI analysis method according to the orcsciit first 
embodiment. FIGS. 5A through 5D are illustrations showing the principle 
imderlying the processing. In a netlist 401, a circuit, which is aH -oblee t the 
subject of the E MI analysis^, is represented as circuit data. Delay information 

20 405 concerning f or e ach node is for m ed den y ed from the netlist 401 through 
using s tatic delay calculation 403 (see FIG. 5A). Transition probabiUty 
information 406 coaccming l or each node is formed derived f rom both the 
netlist 401 and input transition probabiUty 402r -thi-o^igh iismg p ropagation 
probabiUty 404 (see FIG. SB). In consideration of Based on a triangular 
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waveform whose area corresponds to the quantity of electric current derived hy 
means of multiplying current waveform information hy prohability information, 
the average current waveform calculation means 408 derives a n average current 
waveform 409 is formed by average current waveform calculation snean s 408 
5 from element current waveform information 407 eoncemfa Fig -for e ach node (see 
FIG. 5C) and the delay information 405. The thus-determined average current 
waveform 409 is taken u sed a s average current waveform information (see FIG. 
5D). The average current waveform information is subjected to FFT 
processing 410, thereby deriving frequency characteristic information 411. 

10 FIG. 6 shows a flowchart of processing e^- pgrformed by t he average 

current waveform calculation means 408. The average current waveform 
calculation means 408 reads element current waveform information from a 
table (step 1250) and performs a current waveform calculation loop (step 1251). 
The base^ of a triangular waveform of an instance to be processed is 

15 extracted from an output slew (step 1252). The area of the triangular 

waveform is tal^en^ a& beiiig" derived by m<?aHs of multiplying 

transition probability per cycle, and I is taken as the value of the area of the 
triangular waveform. The height."-^%^„_ftj. of the triangular waveform is 

2x1 

calculated from transition probability per cycle (step 1253), wherein-^^ 

20 /^denotes the quantity of electric current flowing in a cell of an event which is an 

object of- i>roce^si^> gbcing processed . This processing corresponds to 

processing performed by a triangular waveform shaping section. Until 
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variftMc ■ x„changes from 0 to W/2, A(c, 0 expressed (by Eq,5) i s added to 
I(t+x) and I(t'X). Further, At is added to variaMe (steps 1254 and 1255). 

Here, /(r+x) denotes total electric current flowing through all the cells at time 
t+Xy and lit'X) denotes total electric current flowing through all the cells at time 
5 t'X. 

The frequency characteristic of a circuit to be Rnalv^ged s ubiecl circuit 
can be determined in the manner as mentioned previously, and a designer can 
analyze EMI wl>i<:4> tliat would arise in a circuit of interest. 

According to the EMI analysis method, a current waveform is modeled 
10 through an averaging operatio n^ on the basis of based on static delay 
information and propagation probabiUty information. The thus-obtained 
model is subjected to FFT processing, thereby analyzing EMI of a circuit. The 
EMI analysis method can analyze EMI components within a shorter period of 
time than e^n-a known gate-level dynamic analysis method cau. 
15 fa a - ease - wber€"perfoiTO@n€e of Wh£B EMI analysis for each path of a 

circuit is desired, static delay information eonccraing to r,each path is given. 

In the i^rescnt fi rst e mbodiment, electric currents of all nodes in a 
€iy€mt"t0- be-amf¥%ed-- subiect cixxuit a re added. However, sf>"4oi^g-"afy"|£_the 
number of nodes whose electric currents are to be added is con tro l led a djusted. 
20 in aceot - daace with as approp^^ magnitude of an electric current or 

the frequency of probability, processing time can be shortened further. 
(Second Embodiment) 

Next;j wiy"be"desc?4b€d an EMI analysis method according to a second 
embodiment of the present inventio n wiM be described . As represented by a 
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fio^v<>l>art In the second embodimeiit. as s hown in the flowchart of F IG. 7, Ifee 
py e s ca t EMI analysis method is character i^ ^ ^ current 
waveform estimation means 708 is nsed i n lieu of the average current waveform 
calculation means 408 employed in of_the first embodiment, and utilisdng 
5 random current waveform information is used i n lieu of the average current 
waveform information. In other respects, the EMI analysis method according 
to the present embodiment is identical in configuration with that described in 
connection with the first embodiment. 

FIG. 7 is a block diagram for dc s eribiBg s howing the overall flow of 

10 processing of the EMI analysis method according to the present s econd 
embodiment. FIGS. 8 A through 8D are illustrations showing the principle 
underlying the processing. In a netlist 701, a circuity which is an object tfie 
subject o f EMI analysis^ is represented as circuit data. Delay information 705 
con ee rni n g for each node is formed derived fr om the netUst 701 through u sing 

15 static delay calculation 703 (see FIG. 8 A). Transition probabiUty information 
706 conccrmng f or e ach node is formed derived fr om the netUst 701 and input 
transition probability 702 , through using propagation probabiUty 704 (see FIG. 
8B). O^ti" the "bj^sIS" of" Based on e lement current waveform information 707 
conccming t or each node (FIG. 8C) and operating frequency information 712, 

20 random waveform estimation means 708 produces random current waveform 
information 709 (see FIG. 8D) within a plurality of predetermined cycles. The 
thus-produced random current waveform information 709 is subjected to FFT 
processing 710, thereby deriving frequency characteristic information 711, 
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FIG. 9 shows a flowchart of processing o ^bv t he random current 
waveform estimation means 708. The average current waveform estimation 
means 708 reads element current waveform information from a table (step 1280) 
and performs a current waveform calculation loop (step 1281). The average 
5 current waveform estimation means 708 performs loop processing until valuable 
^^^» ^the value of v finitiailv 1) changes from 1 to the value of r eaches a giveu 
frequency, (step 1282). The following processing is iterated until calculation of 
a current waveform is completed. A determination is made as to whether m- 
$%et-a random number is smaller than the value of probability (step 1283). If a 
10 random number is smaller, the base of a triangular waveform of an instance to 
be processed is extracted from an output slew (step 1284). At this time, the 

area of the triangular waveform is defined as Wx^, and / is the value of the 
area of the triangular waveform. The height ^%^^of the triangular waveform 

is calculated by 2x— (step 1285), wherein ^^^/denotes the quantity of electric 

W ™ 

15 current flowing in a cell of an events which is an object rhe subject of processing. 
This processing corresponds to processing performed by a triangular waveform 
shaping section. Until variabie ^^^-^ t he value ofx (itiailv 0) c hanges from 0 to 
reaches h(c, i) expresse d ( b y Eq .5) is re peatedly added to I(t+x) and I(t-x). 
Further, AT is added to variable ''x'' x _(steps 1286 and 1287). Here, /(r+x) 

20 denotes total electric current flowing through all the cells at time r+x, and I(t-x) 
denotes total electric current flowing through all the cells at time t-x. 
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The frequency characteristic of a circuit to he axmly/,ed s ubject circuit 
can be determined in the manner a s m co tlo m d previously mentio^^ and a 
designer can analyze EMI wMgfe that would arise in a circuit of interest. 

According to the preseut E MI analysis method^ a current waveform is 
5 modeled through M l^l£„^ random current waveform operatioU j - on the ba^i& - of 
based on s tatic delay information and propagation probability information. 
The thus-obtained model is subjected to FFT processing, thereby analyzing EMI 
of a circuit. The EMI analysis method can analyze EMI components with high 
accuracy within i n a shorter period of time than a known gate-level dynamic 
10 analysis method. 

In a casci where W hen p erformance of EMI analysis for each path of a 
circuit is desired, static delay information eonccrning f or each path is given. 

In the iH^-eseM - sccoml e mbodiment, electric currents information of -for 
all nodes in a circuit to be analya^cd subject circuit are added. However, so long 
15 as-if the number of nodes whose electric currents are to be added is coatrolled 
adjusted, as aoorooriate for i ?> aceordggiee wilih t he magnitude of an electric 
current or the frequency of probability, processing time can be shortened 
further. 
(Third Embodiment) 

20 An EMI analysis method according to a third embodiment of the 

present invention will now be described. In the previously^described first and 
second embodiments, delay information and probability information are 
prepared separately. Information is derived by means of multiplying 
waveform information, which is obtained as element current waveform 
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inforination, by probability information. The thus-obtained information is 
added to a delay time of each node. In contrast, in the present t hird 
embodiment, delay propagation probability information is formed derived f rom 
delay propagation probabiUty information. Delay/transition probability is 
5 calculated from the delay propagation probability information, and element 
waveform information is added to the thus-calculated delay/transition 
probabiUty. 

In this way, more reaUstic current waveform information is calculated^ 
and. ^fee -The r esult of this c urrent waveform calculation is subjected to FFT 

10 processing, thereby determining the frequency characteristic of an EMI 
component of a circuit to be analyzed. Thus, EMI of the circuit is analyzed. 
As can be seen from an enlarged view shown in FIG. 10, the present embodiment 
is directed particularly a case where a pluraUty of paths ar^- is p rovided in a 
composite cell. FIG. 11 shows delay transition information concerning 

15 propagation of a signal in each of the paths of the composite cell shown in FIG. 
10. FIG. 11 shows delay transition probability information as one example. 
As can be seen from FIG. 11, there is obtained n ode informatio n is obtained 
€oncieyni4>g" for a pluraUty of paths, and mean current waveform information is 
for$.Hed derived f rom the node information. 

20 FIG. 12 is a block diagram for describing s ho\^ing the overall flow of 

processing of the EMI analysis method according to the present embodiment. 
FIGS. 13A through 13C are illustrations showing the principle underlying the 
processing. FIG. 14 is a flowchart of average current waveform calculation 
means used in the processing. In a netUst 901, a circuity which is an object t he 
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subject o f EMI analysis* is represented as circuit data. Delay/transition 
probability 906 of each node is calculated from the netlist 901 and input 
transition probability 902----fHV""tl>e--- -1^ on d elav/propagation 

probability 904 (see FIG. 13A). Mean current waveform estimation means 908 
5 produces mean current waveform information 909 (see FIG. 13C), in 
consideratioH of b ased on a triangular waveform whose area is determined by 
the quantity of electric current, such that the delay/transition probability 906 is 
multiplied by element current waveform information 907 (see FIG. 13B). The 
thus- formed calculated m ean current waveform information 909 is subjected to 
10 FFT processing 910 within a time domain wM e h - that i s determined on the b as is 
ef -using o perating frequency information 912, thereby obtaining frequency 
characteristic information 911. 

FIG. 14 shows a flowchart of processing ef -foy t he average current 
15 waveform calculation means. The average current waveform calculation 
means reads element current waveform information from a table (step 1310) 
and performs a current waveform calculation loop (step 1311). The following 
processing is iterated until calculation of a current waveform is completed. 
The delay/transition probability 906 calculated from delay information and 
20 transition probability information is multipUed by element current waveform 
information 907 (see FIG. 13B) (step 1312). Ii> considerati o n of Based on a 
triangular waveform whose area is determined by the quantity of electric 
current, average electric current waveform estimation means 908 adds the 
result of midtiplication as mean current, thereby deriving average current 
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waveform information 909. The average current waveform information 909 is 
subjected to FFT processing 910, thereby determining frequency characteristic 
information 911. 

5 The frequency characteristic of a circuit to be anal yse d sub|ect d^^ 

can be determined in the manner as mentioned previously, and a designer can 
analyze EMI ^ which t hat w ould arise in a circuit of interest. 

According to the present E MI analysis method, delay propagation 
10 probabiUty information is fonn c d derived f rom static delay information and 
propagation probability information, and average current waveform 
information is formed — derived fr om the delay propagation probability 
information. The thus-obtained average current waveform information is 
subjected to FFT processing, thereby enabling highly'-accurate EMI analysis. 
15 The EMI analysis method can analyze EMI components within a shorter period 
of time than ea^a known gate-level dynamic analysis method can . 

In addition to a distribution taking into consideration a path, a 
temperature/process/voltage distribution may be conceived as the 
20 delay/transition probability information shown in FIG. 13 A. 

1*1 each of t he foregoing embodiments, FFT processing has been used for 
analyzing a frequency. However, the present invention is not Umited to FFT 
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processing. Needless to say, another processing method, such as ordinary 
Fourier transformation, may alternatively be employed. 

5 The present invention can mat<?ri8lize - faciMtate e valu of EMI 

developing in an feS¥ -LSI t hrough a simulation, by means o f M^blv ■ very 
accurately analyzing* through static processing, variation in power-supply 
current which may be said to primarily account for EMI. Further, in contrast 
with dynamic analysis of EMI^ which is embodied by a gate-level simulation or a 
10 like simulation, the present EMI analysis method can prevent an increase in 
processing time. 
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